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Abstract 

The motion of a system of particles under electromagnetic interac- 
tion is considered. Under the assumption that the force acting on an 
electric charge is given by the sum of the electromagnetic fields pro- 
duced by any other charged particles in its neighborhood, we prove 
that the vector potential of the electromagnetic field has to be con- 
sidered for the balance of kinetic momentum. The theory cannot be 
quantized in the usual form — because it involves a mass matrix that 
depends on spatial variables — and the Hamilton's function becomes 
singular at a distance equal to the geometric mean of the electrody- 
namic radiuses of electrons and protons. 

1 Introduction 

In previous works (Jj El' have shown that, according to classical 
mechanics and electrodynamics: a neutral system of electric charges 
that passes through a region where there is an inhomogeneous mag- 
netic field, experiences a force, even if its internal kinetic angular mo- 
mentum is equal to zero. Given that this challenges the common in- 
terpretation of the Stern-Gerlach experiment — as evidence that there 
are intrinsic angular momenta — we have considered necessary to study 
the motion of systems of electric charges where the internal magnetic 
force is not neglected, as it is usual in common classical treatments. 
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In this paper we study the motion of electric charges under electro- 
magnetic interaction. Neglecting radiative effects, we assume that the 
force acting on an electric charge is given by the sum of the electro- 
magnetic fields produced by any other particles in its neighborhood. 
Prom the invariance of the Lagrange's function, we find that the vec- 
tor potential of the electromagnetic field must be considered for the 
balance of linear and angular momentum, thus predicting a classical 
Bohm-Aharonov Effect. 

A Hamilton's function is obtained also for the system of two par- 
ticles. The result is a theory that cannot be quantized but approxi- 
mately in the usual form, since it involves a mass matrix that depends 
on spatial variables. Also, the Hamilton's function becomes singular 
where the distance between the particles satisfies the relation: 



For the sake of completeness, we include a section where the La- 
grange's function and the equations of motion for the center of mass 
and the vector of relative position are obtained. 



2 The General Law of Motion 

We study the classical motion of an electron and a proton, under 
electromagnetic interaction. 

Neglecting any retardation and/or radiative effects, we use the 
formulas 



to find the electrodynamic potentials associated to a punctual charge 
q moving along the path f{t). (Where v = The corresponding 
electromagnetic field is 



1 a A 

E(x-t) = -V*--^= (2) 

q{x — r) q I v {{x — r) ■ v)v\ 

Hix,t)=VxA = ^^0^. (3) 
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Further, we suppose — as it's done when only Coulomb's field is 
considered — that the force that acts on the electron is that due to 
the proton's electromagnetic field, and vice versa. The equations of 
motion are: 

m^v, = -'^^^^+'4(^+^-^^^^^^] (4) 

o2 



Ve X {Vp X (fe - fp)) 



and 



Vp X {Ve X {Vp - fe)) 

(Here we have introduced the notation 

r = \\fp - fell (6) 

that simplifies the equations.) 

These are the Euler-Lagrange's Equations for the Lagrange's Func- 
tion: 

T f ^ ^ ^ ^ \ ^ 2 : ^ ^2i^^ ^ Vf.- Vp , . 

L{re,rp,Ve,Vp) = -VfleV^ +-mpV H 2 ^ iV 

2 2 ' r c"^ r 

as we'll prove for the equation 

d dL dL 



dt dvp dfp 



(8) 



From El we get 



d dL _i d f V(. s 



dtdvp~ c'^dtKr 
Further: 

dfVe\_Ve {{rp - fe) ' fp)fe ^ ((^p - ^e) " Ve)Ve ^-j^q^ 



dt \ r J r 
and 

dL _ e'^{rp-fe) ^ e"^ {ve ■ Vp){fp - Ve) 
dfp (? 
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From equations I51ITUI andlTTl — and the identity ax (bx c) = {a-c)b — 
{a - b)c — we can easily prove that equation |H1 is equivalent to equation 

ini 

The generalized momenta are 

€^ V 6 

Pe = meVe = TUeVe Ap{re) (12) 

and 

(? V 6 

Pp = nipVp- = mpVp + -Aeifp) (13) 

where A(>{r) and Ap{r) are the vector potentials for the field of the 
electron and the field of the proton, respectively. 

The Lagrange's function [7| is invariant under translations and ro- 
tations of the reference system; therefore, the sum of the generalized 
momenta 

€ -* € -* 

Pe+Pp = rUeVe + rUpVp Ap{fe) + -^e(rp), (14) 

c c 

and the total angular momentum (which is not equal to the kinetic 
momentum and, therefore, the dipolar field is not enough to describe 
the magnetic properties of the system) 

L = reXpe + rpX pp, (15) 

are constants of motion. In consequence, the center of mass of the 
system, 

j^^ m,r, + mpfp ^ (16) 
rUe + mp 

does not move according to Newton's First Law. This was expected 
given that equations 0& El are not in compliance with Newton's Third 
Law either. 

In the case of a system of n particles with masses mi, • ■ ■ , and 
charges qi,- " ^Qn, the Lagrange's Function assumes the form 

i = i:r'"t^-iEf^(i-^). a^) 

i=l ^ ^ q,^q, ^ ^ ^ 

Therefore, the generalized momentum of the ith particle is 



Pi = niiVi H Y] — — = rrnvi H Ai{ri), (18) 

c c n,- 
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where Ai is the vector potential of the magnetic field produced by the 
other particles. 

Given that the function 1171 is also invariant under arbitrary trans- 
lations and/or rotations, we come again to the conclusion that the 
sum of the generalized momenta and the angular momentum, are con- 
stants of motion. Therefore, as it has been confirmed by Bohm and 
Aharonov 3j, the vector potential of the electromagnetic field acting on 
each particle must be considered for the balance of linear momentum. 

The problem of gauge invariance is not an issue for us. Under the 
gauge transformation 

the momenta are transformed as: 

Pi' = p, + ^, (20) 

If the sum of the momenta is going to be a constant of motion, A must 
be invariant under arbitrary translations. In other words 



pi) 



• 1 dr, 
1=1 



and, in those circumstances: 

n 



n 



i=l 1=1 

Going back to the electron-proton system, its energy 

E = ^ruevj" + -m,^/ - - - (22) 

1,. . . e2 

= 7;{Pe -Ve+Pp- Vp) , 

1 r 
is also a constant of motion. 

Solving [T^ and IT^ for the velocities, we find 



From this and 1221 we get the Hamilton's Function 



Pe^ _|_ _ e^VePp 2 



^ = ^ /"-'"^^ ' _ (25) 



Where 



e2 



^ » ^7 u72 (26) 

this function coincides with the function used in W to approximate 
the eigenvalues of the corresponding quantum system. 
Further, we notice that ESI is singular where 



e2 



c^(memp)-'^/^ 



(27) 



At this distance the equations of motion ^ and IS] cannot be solved 
for the accelerations, which is fundamental for the applicability of the 
theorem of existence and uniqueness of solutions. Therefore, even 
if the Principle of Least Action is valid, to determine a particular 
solution additional conditions have to be imposed. There are another 
two possibilities which we shall not investigate further in this paper: 

1. That inequality defines the limits of validity of electrodynam- 
ics. 

2. That a fully relativistic approach is required. In this case the 
effects of retardation have to be considered; the equations of mo- 
tion are difference-differential equations; and there is not room 
for a variational approach. (At least not for a variational ap- 
proach that does not explicitly accounts for the action of the 
entire electromagnetic field.) 



3 Separation of the Internal Motion 

As it was shown before, the center of mass of the system does not 
move according to Newton's First Law. Notwithstanding, and for the 
sake of completeness, we'll carry out the decomposition of the motion 
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into the motion of the center of mass and an internal motion. Let's 
consider the substitutions: 

R=—tLJL— ; r = re-rp (28) 

(where M = nip + nie), in such way that: 

rp = R- -^r- fe = R + —^r. 29 



Vp ■ Ve = R + KlR ■ f Ty^o ^ 



and 

5^ I 5 - '^errip-^ 

where 

nip - nie 

The Lagrange's function [7| takes the form: 

R +1 

L{R,f,R,r} = ^MR +^^f=" + — -^-I M±_, (30) 



2 ' 2""' ' r 



The momenta are: 



and 

KLe"" A 2e2 



^>=-^^ + Hi + ma: - (32) 



The energy is: 



.2 ^ . ^ 



e2 e^i? + KLR-r-^^ r 
r 

„2 



E{R,f,R,f) = -MR +-/xr ^ Ml (33) 

2 2 r r 



= ^(P^.i?+p>-rl--. 
Solving equations 1^ and for the velocities we find 



and 



A 



(35) 
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where 

A = M/u I 1 



\ Kle^ _ 



Now we are ready to write the Hamilton's Function 



2 ^ e4 7 (36) 



rripmeC^r^ 
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